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The large disparity between the scale of quantum grav- 
ity, i.e., the Planck scale, Mp\ » 10 16 TeV, and the elec- 
troweak scale, of the order of 1 TeV, is known in the 
standard model (SM) as the hierarchy problem. In the 
presence of this hierarchy of scales it is not possible to 
stabilize the Higgs boson mass at the low values required 
by experimental data, unless by using an unlikely large 
amount of fine-tuning. 

In the Randall-Sundrum model [1], the existence of a 
fifth dimension with a warped spacetime metric is pro- 
posed, bounded by two three-dimensional branes. The 
SM fields are localized on one brane, while gravity orig- 



inates on the other. With this configuration, TeV scales 
are naturally generated from the Planck scale due to a 
geometrical exponential factor (the "warp factor"), = 
Mpiexp(— k-Kr c ), if kr c m 12, where Mpi = Mpi/y/sn is 
the reduced Planck scale, and k and r c are the curvature 
scale and compactification radius of the extra dimension, 
respectively. 

Gravitons are the only particles that propagate in the 
fifth dimension, and appear as a Kaluza-Klein series 0] of 
massive excitations (KK gravitons, G) with spin 2, mass 
splitting of the order of 1 TeV, and a universal coupling 
to the SM fields. Phenomenologically, it is convenient 
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to express the two Randall-Sundrum parameters k and 
r c in terms of two direct observables: the mass of the 
lightest excitation, Mi, and the dimensionless coupling 
to the SM fields, fc/M P1 . To address the hierarchy prob- 
lem without the need for fine-tuning, Mi should be in 
the TeV range and 0.01 < fc/M P1 < 0.1 [3]. KK gravi- 
ton resonances could be produced in high energy particle 
collisions and would subsequently decay to pairs of SM 
fermions or bosons. 

In this Letter, we report an inclusive search for the 
lightest KK graviton in the ee and 77 decay channels 
with the DO detector at the Fermilab Tevatron Col- 
lider, where protons and antiprotons collide at ^/s = 
1.96 TeV. KK gravitons would be produced via quark- 
antiquark annihilation and gluon-gluon fusion processes. 
For fc/Mpi < 0.1, KK gravitons would appear as nar- 
row resonances in the ee and 77 invariant mass spectra, 
with a natural width much smaller than the resolution 
of the DO detector and with a branching fraction for the 
77 decay mode which is twice that of the decay to ee. 
Previous DO searches for KK gravitons have excluded 
Mi < 300 GeV for fc/M P1 = 0.01 and Mi < 900 GeV 
for k/Mpi = 0.1 at the 95% C.L. A]. CDF has recently 
excluded Mi < 889 GeV for fc/M P1 = 0.1 at the 95% 
C.L. d. 

The DO detector @, 0] consists of tracking detectors, 
calorimeters, and a muon spectrometer. The tracking 
system includes a silicon microstrip tracker close to the 
beam and a central fiber tracker, both located within 
a 2 T superconducting solenoidal magnet. The liquid- 
argon and uranium calorimeters consist of a central sec- 
tion covering pseudorapiditics \r]\ < 1.1 and two end 
cap calorimeters that extend the coverage to \rj\ s» 4.2, 
where 77 = — ln[tan(6>/2)], and 9 is the polar angle with 
respect to the proton beam direction. The azimuthal 
angle is denoted by 4>. The electromagnetic (EM) sec- 
tion of the calorimeters is segmented into four longitu- 
dinal layers (EMi, «=1,4) with transverse segmentation 
of Arj x A<fi = 0.1 x 0.1, except for the more finely seg- 
mented EM3 section where it is 0.05 x 0.05. A preshower 
system (CPS) uses plastic scintillators with different ori- 
entations located between the solenoid and the cryostat 
of the central calorimeter and provides precise measure- 
ments of the positions of EM showers. The luminosity is 
measured using plastic scintillator arrays placed in front 
of the end cap calorimeters. The data sample was col- 
lected between July 2002 and June 2009 using triggers 
requiring at least two clusters of energy deposits in the 
EM calorimeter and corresponds to an integrated lumi- 
nosity of 5.4 ±0.3 fb -1 . 

We select events with two EM clusters, each with trans- 
verse momentum pr > 25 GeV and \r]\ < 1.1, recon- 
structed in a cone of radius R = yj (A-q) 2 + (Acj)) 2 =0.4. 
The EM clusters are required to have at least 97% of 
their energy deposited in the EM calorimeter and to 
have the calorimeter isolation variable / = [Et o t(0A) — 



^em(0.2)]/^em(0.2) < 0.07, where E tot (R) [E EM (R)} is 
the total [EM] energy in a cone of radius R. 

Given the different branching fractions for the 77 and 
ee decays of the KK graviton, plus the fact that the two 
channels have different backgrounds, the analysis treats 
the two channels separately to optimize the sensitivity. 
If both EM clusters in an event are spatially matched to 
tracker activity, either a reconstructed track or a den- 
sity of hits in the silicon microstrip tracker and central 
fiber tracker consistent with that of an electron, the event 
goes in the ee category. Otherwise, the event is put in 
the 77 category, which contains events with at least one 
EM cluster failing to match tracker activity. With this 
definition, about 97% of the selected G — > ee events are 
put in the ee category and about 90% of the selected 
G — » 77 events are put in the 77 category. 

In the ee category, the two electrons are not required 
to have opposite charges to avoid the loss due to charge 
misidentification, and two additional requirements are 
placed on each EM cluster: (i) the scalar sum of the 
Px of all tracks originating from the primary vertex (PV, 
see below) in an annulus of 0.05 < R < 0.4 around the 
cluster, I tr k, be less than 2.5 GeV; (ii) the cluster be con- 
sistent with the electron shower shape using a x 2 test an d 
a neural network discriminant Q . In the 77 category, ad- 
ditional requirements are placed on each EM cluster: (i) 
Itrk < 2.0 GeV; (ii) the energy-weighted shower width 
in the r — <fi plane in EM3 be less than 3.7 cm; (hi) the 
cluster be consistent with the photon shower shape using 
a neural network discriminant. 

Proper reconstruction of the event kinematics requires 
correct identification of the PV of the hard collision. For 
events in the ee category, the PV is chosen from the list 
of vertex candidates as the one with the least probability 
of being a vertex from a soft pp interaction as estimated 
from the pr of associated tracks. For the 77 category, we 
use the EM-CPS pointing capability, which reconstructs 
the axes of EM showers by fitting straight lines to shower 
positions measured in the four longitudinal calorimeter 
layers and the CPS. The EM-CPS pointing spatial reso- 
lution is 3.7 ±0.2 cm along the beam axis. If at least one 
photon candidate is matched to a CPS cluster the 
vertex consistent with the EM-CPS pointing position is 
chosen as the PV. For events with no photon candidate 
having a CPS match or events with inconsistent EM-CPS 
pointing positions of the two photon candidates, the PV 
is chosen as the one with the highest number of associ- 
ated tracks. The PV is required to lie within 60 cm of the 
geometrical center of the detector along the beam axis. 
The data include a total of 203586 events (186596 in the 
ee category and 16990 in the 77 category) that satisfy 
these selection criteria and with the invariant mass of 
the two EM clusters M ee / 77 > 60 GeV. 

All Monte Carlo samples used in this analysis were 
generated using pythia [10| with CTEQ6L1 [ll[ parton 
distribution functions, and processed through a GEANT- 



5 



based [12] simulation of the DO detector and the same 
reconstruction software as the data. KK graviton sig- 
nals in the ee and 77 decay channels are simulated over 
the range of parameters 220 < M\ < 1050 GeV and 
0.01 < fc/Mp] < 0.1. The accuracy of the PV asso- 
ciation has been studied in KK graviton events, where 
the PV reconstruction efficiency is ss 98%, with w 96% 
(« 93%) probability to match the true vertex in the 
ee (77) channel. The simulated and observed invariant 
mass spectra are compared in ee and 77 categories sep- 
arately. The dominant irreducible background in the ee 
final state is due to the Drell-Yan (DY) process, where 
an ee mass-dependent k factor [l3l ] is applied to correct 
the pythia spectrum for next-to-next-to-leading order 
effects. The dominant irreducible background in the 77 
final state is SM 77 production, where pythia 77 events 
are reweighted to reproduce the 77 invariant mass spec- 
trum predicted by the next-to-leading-order calculation 
of diphox QJJ. DO has measured the SM 77 differen- 



tial cross section with respect to the 77 invariant mass, 
and in the range used for this analysis (above 60 GeV) 
the shape of this distribution from diphox agrees with 
the data 15] . The leading systematic uncertainty on this 



background's shape arises from the choices in the scales 
used in the diphox calculation, and is at the level of 
10%. The main instrumental background comes from 
the misidentification of one or two jets as electrons or 
photons. The shape of the invariant mass spectrum of 
this source of events is estimated from data by selecting 
events with EM clusters that are not consistent with elec- 
tron or photon showers using the y 2 test (ee category) 
or the neural network discriminant (77 category). Other 
SM backgrounds, due to DY rr, IF + 7, WW, ZZ, WZ, 
VF+jets, and ti production, are small and are estimated 
using pythia Monte Carlo events corrected to account 
for higher order effects lB- 18j|. 



Having obtained the shapes of the invariant mass spec- 
tra of the various background sources, the background 
normalization is determined by fitting the invariant mass 
spectrum of the data to a superposition of the back- 
grounds in a low-mass control region (60 < M ee / 77 < 
200 GeV) , where KK gravitons have been excluded at the 
95% C.L. by previous searches. In the fit, the total num- 
ber of background events is fixed to the number of events 
observed in the data, and the contributions from SM 77, 
DY ee, and instrumental background are free parame- 
ters, while the other SM backgrounds are normalized to 
their theoretical cross sections. The fit is performed for 
the ee and 77 categories separately. By varying the crite- 
ria to select the instrumental background sample and the 
fitting range, the uncertainty of the background normal- 
ization procedure is estimated at the level of 2% (10%) 
in the ee (77) category. 

Figure [T] shows the measured ee and 77 invariant 
mass spectra from the data, superimposed on the ex- 
pected backgrounds. The data and predicted background 
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FIG. 1: Invariant mass spectrum from (a) ee and (b) 77 
data (points). Superimposed are the fitted total background 
shape from SM processes including instrumental background 
(open histogram) and the fitted contribution from events with 
misidentified clusters alone (shaded histogram). The open 
histogram with dashed line shows the signals expected from 
KK gravitons with Mi = 300, 450, 600 GeV (from left to 
right) and fe/Mpi = 0.02 on top of the total background. In- 
variant masses below 200 GeV are taken as the control region. 



are generally in good agreement. In the region around 
450 GeV there is an excess of events in the 77 invari- 
ant mass spectrum. As estimated with pseudoexper- 
iments, the probability that this excess is exclusively 
due to backgrounds' fluctuations is 0.011, implying that 
the background-only hypothesis is disfavored at the 2.30 
standard deviations (s.d.) level. If we assume that this 
excess is due to a KK graviton, including the ee channel 
reduces the significance to 2.16 s.d.. 

In the absence of any significant signal for a heavy 
narrow resonance, we compute upper limits for the pro- 
duction cross section of KK gravitons times the branch- 
ing fraction into the ee final state using a Poisson log- 
likelihood ratio (LLR) test [l9j]. Invariant mass distri- 
butions are utilized in the limit calculation. The ee and 
77 categories are treated as two independent channels, 
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and then the two separate LLRs are added to obtain a 
combined exclusion limit assuming the 1:2 ratio of the 
branching fractions. 

Systematic uncertainties on the backgrounds' predic- 
tions and on the signal efficiency are considered to cal- 
culate limits. These include the integrated luminosity 
(6.1%), parton distribution functions (0.7% - 6.6% for 
the acceptance and 9.2% - 16.9% for the graviton pro- 
duction cross section), electron and photon identification 
efficiency (3.0% per object), EM cluster energy resolution 
(6%), and trigger efficiency (0.1%). The uncertainty on 
the acceptance due to initial state radiation (ISR) is esti- 
mated to be 4% by varying the parameters governing ISR 
in pythia. Uncertainties affecting the expected back- 
grounds arise from electron and photon identification ef- 
ficiency (3.0% per object), mass dependence of the DY 
ee next-to-next-to-leading order k factor (5.0%), shape 
of the SM 77 invariant mass spectrum, and background 
normalization. For the EM energy resolution, the SM 
77 invariant mass spectrum and the background normal- 
ization we consider both the effects on the normaliza- 
tion and on the shape of the invariant mass distribution 
used in extracting limits. For all other systematic sources 
we consider only changes to the overall background nor- 
malization or signal detection efficiency. Systematic un- 
certainties are incorporated via convolution of the Pois- 
son probability distributions for signal and background 
with Gaussian distributions corresponding to the differ- 
ent sources of systematic uncertainty. Correlations in the 
systematic uncertainties between signal and background 
in ee and 77 categories are taken into account. 

The resulting limits on the production cross section 
times branching fraction into electron-positron pairs of 
the lightest KK graviton, cr(pp —> G + X) x B(G — > ee), 
are given in Table Q] and displayed in Fig. [2] As shown 
in Fig. [3l using the cross section predictions from the 
Randall-Sundrum model with a k factor of 1.54 J2fj, we 
can express these results as upper limits on the coupling 
k I Mpi as a function of M\ . 

In summary, using 5.4 fb -1 of integrated luminosity 
collected with the DO detector at the Fermilab Teva- 
tron Collider, we have searched for a heavy narrow res- 
onance in the ee and 77 invariant mass spectra. The 
observed spectra agree with predictions from SM back- 
ground processes. For the Randall-Sundrum model with 
a warped extra dimension, we set 95% C.L. upper limits 
on a(pp — >• G + X) x B(G —> ee) of the lightest Kaluza- 
Klein mode of the graviton between 6.7 fb and 0.43 fb 
for masses between 220 and 1050 GeV at the 95% C.L., 
which translate into lower limits on the mass Mi of the 
lightest Kaluza-Klein excitation of the graviton between 
560 and 1050 GeV for couplings of the graviton to the 
SM fields 0.01 < fc/Mpi < 0.1. These results represent 
the most sensitive limits to date. 
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FIG. 2: 95% C.L. upper limit on a(pp -> G + X) x B(G -»• 
ee) from 5.4 fb - of integrated luminosity compared with the 
expected limit and the theoretical predictions for different 
couplings fc/Mpi. 



TABLE I: 95% C.L. upper limit on a(pp -> G + X) x B(G -> 
ee) and coupling k/Mpi from 5.4 fb -1 of integrated luminos- 
ity. 



Graviton Mass 


a x B(G - 


-> ee) (fb) 


Coupling fc/Mpi 


GeV 


Expected 


Observed 


Expected 


Observed 


220 


10.62 


6.71 


0.0034 


0.0027 


250 


7.18 


5.23 


0.0038 


0.0033 


270 


5.91 


5.69 


0.0042 


0.0041 


300 


4.00 


5.37 


0.0044 


0.0050 


350 


2.67 


3.30 


0.0051 


0.0056 


400 


2.12 


1.52 


0.0062 


0.0053 


450 


1.40 


3.03 


0.0068 


0.0099 


500 


1.15 


1.31 


0.0081 


0.0087 


550 


0.89 


0.90 


0.0093 


0.0094 


600 


0.75 


0.84 


0.0111 


0.0117 


650 


0.65 


0.68 


0.0133 


0.0136 


700 


0.56 


0.48 


0.0160 


0.0147 


750 


0.53 


0.52 


0.0199 


0.0197 


800 


0.48 


0.48 


0.0248 


0.0247 


850 


0.46 


0.44 


0.0316 


0.0312 


900 


0.44 


0.43 


0.0406 


0.0403 


950 


0.44 


0.43 


0.0545 


0.0539 


1000 


0.43 


0.43 


0.0713 


0.0713 


1050 


0.43 


0.43 


0.0969 


0.0964 
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FIG. 3: 95% C.L. upper limit on fc/Mpi versus the graviton 
mass Mi from 5.4 fb _1 of integrated luminosity compared 
with the expected limit and the previously published exclu- 
sion contour [J. 
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